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Abstract 

We consider a (long-lived) sterile neutrino dark matter scenario in a five dimensional (5D) 
warped extra dimension model where the fields can live in the bulk, which is partly motivated 
from the absence of the absolutely stable particles in a simple Randall-Sundrum model. The 
dominant production of the sterile neutrino can come from the decay of the radion (the scalar 
field representing the brane separation) around the electroweak scale. The suppressions of 
the 4D parameters due to the warp factor and the small wave function overlaps in the extra 
dimension help alleviate the exceeding fine-tunings typical for a sterile neutrino dark matter 
scenario in a 4D setup. 

PACS: 98.80.Cq 

1 Introduction 

We will discuss the cosmological aspects of the radion phenomenology in the Randall-Sundrum 
background with the fields in the bulk [HElEllllElEllHlEl [ini HI]. In particular we study the decay 
of the radion into the sterile right-handed neutrino which can account for all the dark matter in the 
Universe. This is partly motivated from the absence of the discrete symmetry in a simple warped 
extra dimension model and consequently the lack of an absolutely stable particle which could be a 
promising dark matter candidate. This is in contrast to, for instance, a flat extra dimension model 
which possesses the discrete symmetry corresponding to the translational invariance in the extra 
dimension, called the Kaluza-Klein (KK) parity, leading to the absolutely stable KK dark matter 
II21[I3]0. We consider a simple addition of three right-handed neutrinos to the minimal Standard 
Model particle contents, where one of them is sterile (i.e. very weakly coupling to the other fields) 



^Ref. '14] investigates the dark matter in a warped extra dimension model, where the discrete symmetry was 
introduced leading to the absolutely stable dark matter candidate. 
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with a lifetime longer than the age of the Universe (and the other two explain the atmospheric and 
solar neutrino data [15] for the finite small neutrino masses). 

The sterile neutrino can be either warm or cold dark matter in our scenarios. 

The sterile neutrino warm dark matter in such a minimal extension of the Standard Model 
(sometimes referred to as i/MSM [161 [13 [13 [I9] in a 4D setup) would have a potential interest 
from the astrophysical viewpoint because it can ameliorate the shortcomings in the small (i.e. 
galactic) scale structures of the cold dark matter scenarios which the N-body numerical simulations 
apparently suffer from, such as the missing satellite problem (the cold dark matter models predict 
too many satellites (small dwarf galaxies about a thousandth the mass of the Milky Way) than 
observed) and the cusp/core problem (the cusped central density distributions rather than the 
observed smoother core in the dark matter haloes) [20] . 

One of the notable features in a warped extra dimension setup is the production mechanism of 
the sterile neutrinos which can come from the decay of the radion because the radion can couple to all 
the degrees of freedom. This would be of significant interest for a warm dark matter scenario where 
a simple Dodelson-Widrow ((non- resonant) active-sterile neutrino mixing [2Tl[22]) mechanism alone 
cannot account for all the dark matter of the Universe due to the conflicts from the astrophysical 
observations such as the Lyman-a forest (giving the sterile neutrino dark matter mass lower limit 
of at least ~ 10 keV) and the X-ray data (giving the upper bound of at most ~ 8 keV) [23l [211 [SSI 



The big fine-tunings, such as those typical for a sterile neutrino dark matter scenario in the 4D 
setups, are relaxed in the warped extra dimension models because the 4D parameters are exponen- 
tially suppressed naturally due to the warp factor and the small wave function overlaps in the extra 
dimension. We also mention that the desirable baryon asymmetry of the Universe in our simple 
warped extra dimension model can arise from the decay of the higher KK modes of the right-handed 
neutrinos [52] . 



We consider the scenario where the fields can live in the bulk of the 5D Randall- Sundrum spacetime 
whose AdS^ metric in the conformally flat coordinate is [T] 



[261 [271 EH Ea [301 [31]. 
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where z is in the interval [-R, B!\ with the UV (Planck) brane and the IR (TeV) brane respectively 
at z = R, R' . The radion corresponds to the scalar perturbation F of the metric [5l [TTl [331 El] H 

ds" = (^^^ {e-^^7]^,dx>'dx'' - (1 + 2Ffdz'). (3) 

Assuming the form F{x,z) = f{z)r{x), the canonically normalized radion r{x) in 4D is obtained 
by 

The linear variation of the action 

I d'x^FiTiT'''' - 3T''g,,) (5) 

can give us a coupling of the radion to the matter fields, where g is the determinant of the 5D 
metric and we use the upper-case Latin (lower-case Greek) letters to denote the 5D (4D) indices. 

Let us now consider a concrete example of a 5D bulk Dirac spinor consisting of two two- 
component spinors = (Xo?^") whose standard 5D bulk action including the bulk Dirac mass 
term reads P El [101 [11] 



d^x 

z 



(6) 



where (R/z)^ comes from the determinant and an additional factor (z/R) comes from the 5D 
vielbein . We shall use the standard parameterization of the bulk mass parameter rriu = c/R. 
Note that we did not perform the integration by parts for the 5D kinetic terms because, as we shall 
see in the following section, there arise the non-trivial boundary conditions due to the finite interval 
in the 5th dimension. By using the perturbations in the metric 

V9= (f)'e-^^(l + 2F), ef = diag|(e^e^e^e^l/(l + 2F)) (7) 



^The general scalar perturbations about the background with the vacuum expectation value (vev) of Goldberger- 
Wise (GW) field 0o [11 |3j is (bulk metric scalar fluctuations can be sourced by the perturbation of the GW field as 
well) 

(f>ix, z) = 0o(z) + ^{x, z), ds" = (^^^ (e-^^f^'^^jy^.dx^dx" - (f + G(x, z)fdz^) (2) 

Plugging this parameterization ansatz into the Einstein equations reveals G and (p can be expressed in terms of F 
(in particular G = 2F up to the linear order). Therefore we shall use only a single scalar mode F to describe the 
scalar perturbations in our discussions. 
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the linear order couplings of the radion to a fermion for the above action can be obtained from |f| 



z 



We are mainly interested in the radion couplings to the lightest 4D modes in the KK eigenstate 
decompositions 



X 



^9n{z)Xn{x),-^ = ^fn{z)lpn{x) (9) 



where the 4D components satisfy the 4D equations of motion with the mass m„ for each n^^ mode 

- ia^df^Xn + rrini^n = 0, -ia'^d^tpn + rrinXn = (10) 

In particular, for n = 0, there are massless zero modes of form, with the standard Dirichlet boundary 
conditions 4'\z=r,r' = 0, 



. , , /I - 2c 1 R 

where the normalization constants are obtained by requiring the canonically normalized 4D kinetic 
terms. Eq. (11 II) shows the fermion zero modes are localized around the Planck (TeV) brane for 
the bulk mass parameter c > 1/2 (c < 1/2). The radion, on the other hand, is always localized 
to the TeV brane as can be seen by Eq (jl]). This means that, from the AdS/CFT correspondence 
|36j , the radion is a composite state which shows up when the conformal symmetry is broken at the 
temperature around the TeV scale. We implicitly assume the confinement of the composites when 
we discuss the radion decay around the electroweak scale temperature. 



3 Sterile neutrino mass 

To obtain the abundance of the sterile neutrinos produced by the decay of the radion, we need to 
know the sterile neutrino mass in addition to the couping strength to the radion outlined in the 
last section. 

We assume, for simplicity, three gauge-singlet right-handed neutrinos in addition to the minimal 
Standard Model particles, and consider the possibility for one of them, denoted as A^, to be the 
sterile neutrino dark matter taking account of all the dark matter in the Universe. For it to be 
sterile, we assume the 4D Dirac Yukawa coupling of to the Higgs and left-handed neutrino is 
negligibly small (consequently its Dirac mass is also negligible compared with its Majorana mass) 
which could be justified by the suppression due to the warp factor and small wave function overlaps 

■^We here consider only the hnear order couphng for simphcity. The higher order couphngs such as 
Xn{x)<y'^Xnix)di,r(x) coming from the canonical normahzation of the fermion field [6] and the effects of the ra- 
dion Higgs mixing through the curvature-Higgs couplings |341 I35j would not affect the qualitative analysis even 
though they could potentially affect the quantitative discussions for some parameter range. 
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in the extra dimension (more quantitative discussions on how small it needs to be will be given 
in Because no symmetry prohibits the Majorana mass terms for the gauge singlet A^, we can 
consider the lightest eigenmass of N in existence of the brane localized Majorana mass term which 
dominates over the Dirac mass contributions. 

We briefly review here, for the illustration purpose, the Majorana mass term confined on the 
Planck brane (we choose the basis where the Majorana mass is real) [HI [9l [TOl [STJ [38] 

j d^x^^niMiNRNR + h.c), tum = dwi^iz - R) (12) 

and the boundary conditions are accordingly modified to 

Nl{z = R) = dMNniz = R), Nl{z = R') = (13) 
where dM is a dimensionless constant and 

N{x,z) = ( j^j^'^^j ) ,iV, = $:/r(z)iVf (a:), Nn = fn\')NR\^) (14) 

^ ^ n n 

We can however simplify our analysis by calculating the 4D mass eigenvalues/states via the di- 
agonalization of the mass matrix in the basis obtained without the boundary Majorana mass 
(A^^^-*, nI^\ ...). In this approximation, the symbols /^"|, N^]^ are used for the states obtained 
without the Majorana mass terms rather than the mass eigenstates, and the lightest eigenstate 
dominantly consists of N^^ while A'^^^'' is decoupled from the low-energy theory as the result of the 
boundary conditions. The approximate mass spectrum up to the n^^ KK level can be calculated by 
truncating at the n^^ level mass matrix when the mixing contributions from the higher modes are 
small [371 139] • We are interested in the lightest 4D eigenstate (because of the boundary Majorana 
mass, the lightest mode now obtains a non- vanishing eigenmass ttin), and we omit the upper mode 
indices in the following unless stated otherwise. For instance, the lightest 4D eigenmass for the 
range of our interest m^v -C l/z can be obtained, up to the leading order, by substituting the wave 
functions of Eq. (fTTl) into the brane localized Majorana mass term of Eq. (fT2l) as (note our use of 
c is such that c < 1/2 (c > 1/2) for the localization to the TeV (Planck) brane) [fl 

1 , rR-'-'"^ 



rriM ~ -dM ( ) (1 - 2c) for c < 1/2 (15) 



''The Majorana spinors cannot exist in 5D because of the lack of the real representation of 7 matrices [5, but 
we can add the gauge- and Lorentz-invariant bilinear terms for 5D Dirac isosinglet neutrino, N'^C^^'>~^N (C^^^ is 
a charge conjugation operator), which can be interpreted as the Majorana mass terms in 4D. Note each 5D Dirac 
spinor in general leads to two 4D Majorana states at each KK level in absence of the boundary conditions. 

^ We consider only c < 1/2 because mjv ~ (iM(2c— l)/i?for c > 1/2, and |mArz| <c 1 would require the fine-tunings 
in this case. Note we can obtain the same results as Eqs. ([T5| [T9\\ by substituting the eigenstate solutions (in terms 
of the Bessel functions) of the 5D bulk equations of motion into Eq. (fTS]) and expansing the Bessel functions to the 
leading order for Im^v^:! <C 1. 
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The leading-order radion decay channel into the zero mode sterile neutrinos (approximated by 
N^\x)) can be thus read off from Eq. ([8]) as[§ 



d'xao^rix) < (x)iVf (x) + = j dz (-|) [^j W/f (^) (18) 

The analogous arguments can be applied to the Majorana mass confined on the TeV brane, 
where the lightest eigenmass becomes 



1 , /R^'^-' 



ruN ^ -^dM ( ) (2c - 1) for c > 1/2 (19) 



4 Sterile neutrino abundance 

For an estimation of the abundance of the sterile neutrino dark matter, we integrate the Boltzmann 
equation {Ccoi is a collision term) 

^ + ^Hnj, = (20) 
which can be rewritten in terms of the so-called yield parameter Y = n/s 

dY -C^oi ^ T dgs{T) \ ^21) 



dT HTs V 3c/, (T) dT 

In the above, is the number density for the sterile neutrino, s is the entropy density and H is 
the Hubble constant given by, in terms of the radiation energy p/j, 

H = -^PR^ PR = ^9*T , s = -j^9sT (22) 



^We here note that the radion couples only to the bulk terms and the brane-localized Majorana mass terms do 
not affect the radion couplings [5 . At first sight, one may naively expect that the Majorana mass terms ol Eq. (fT2|) 
could lead to the radion coupling of 

J d^xe^^^^dMiNRNR + h.c.^^R (16) 

This radion coupling term however is exactly canceled by the surface term of the action (see Eqs. ([51 [5])) {[x]^ = 
x{R') - x{R)) 



^ Sboundary — 2 ^ ^ " 



R 
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-I R' 



e-^^{SNLNR - SNrNl + h.c.) 



(17) 



J R 



because of the boundary conditions Eq. p^ . 

^We consider only c > 1/2 for the TeV brane confined Majorana mass case because tojv « dA/(l — 2c) /R' for 
c < 1/2, and |mArz| <^ 1 would require the fine-tunings. 



6 



The total numbers of the effectively massless degrees of freedom are given by 



where gB{F) is the number of helicity states for each boson (fermion) with their corresponding 
temperature Tb(f)- In deriving Eq. lpTI) . we assumed the constant entropy and used Eq. ( |22|) . 

We note that the thermalization of the radion can be justified from the interactions involving 
the gauge couplings. For instance, its coupling to the gauge fields in 4D |H |5] can be obtained in 
the same manner as those outlined in [|2] 

where F is the field strength of the massless gauge field A (photon or gluon) and g^, are the 4D and 
5D gauge couplings which are related via g"^ ~ glR\og{R' / R) . Note that the coupling constant has a 
suppression of the TeV scale which is the cut-off scale on the TeV brane. For example, the interaction 
rate of rA ^ ff around the TeV scale temperature is F ~ (e/ \og{R' /R)Y x TeV (e is the electric 
charge of a fermion /) which indeed exceeds the Hubble expansion rate H ~ {TeV/Mp) x TeV 
(from the 4D Friedmann equation), so that we can use the thermal abundance for the radion. 

Let us now parameterize the collision term as Ccoi = rir'jX'^mr/STT. Here n^, rrir denote the radion 
abundance and mass and A represents a dimensionless (4D) coupling constant between the radion 
and the sterile neutrino (in particular, this parameterization of the collision term corresponds to 
the decay channel XrNN to be discussed as our concrete examples in the next section). 7 is the 
mean time dilation factor in the decay and we use the Maxwell-Boltzmann statistics for the phase 
space distribution function f{p) 

g.T^ fmrY fmr\ ImA J d^Pfjp)^ K,{m.,/T) 



2ti^ \ T J "\TJ' ' \Er/ fd^pfip) K2{mr/T) 

where Kn is the modified Bessel function of the second kind of order n (as expected, rir oc 
g^{mrTy^'^exp{—mr/T),'y ~ 1 for T ^ rrir and Ur oc g^T^,'j ~ m/T for T ^ rur) [lOl 1^. 
Then Y can be obtained by integrating the Boltzmann equation 



nir 



X 



Y{x) = aX^^y{x), y{x) = I dx' Ki{x')x'\ (26) 

Jnir/HTeV 



where x = rrir/T and the numerical coefficient a ~ 2.6 x 10^'^ with g^ ~ 110. The numerical 
integration of Ki{x)x^ approaches quickly to the asymptotic value of y{x ^ 1) ~ 4.7 where we 
took the integration lower bound of rrir/ fiTeV = 0.1 {f^TeV is the scale of the conformal symmetry 
breaking) |^. 



^There is an uncertainty for when the conformal symmetry breaking occurs, even though changing the lower 
bound to rrirl ^TeV = 0.01 changes the asymptotic value less than 1% and m^/ fJ-TeV — 1-0 just slightly changes the 
asymptotic value to y{x » 1) ~ 4.5. 
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In the above integration, the time dependence of gs was assumed to be small. This assumption 
would be reasonable during the decay of the radion which dominantly occurs around T ~ (the 
electroweak scale). The most of the higher KK states which could possibly have a non-negligible 
effect on dgs{T)/dT are not energetically accessible around the conformal symmetry breaking energy 
scale when the composite particles appear. More detailed studies of the physics around the phase 
transition period will be presented in our future work |^. Matching m^Y to the current matter 
density value Pm/ s ~ 0.4 x 10~^GeV |43| gives an estimation of 



A 



2 



0.3 X 10-- ( ( ) (27) 



') \imGev) 



5 Examples 

Let us now discuss a few concrete examples for the illustration purpose. For definiteness, we set 
to be the Planck scale and use the values for the radion mass and Majorana mass parameter 
as ttIt = 300 GeV and (Im = 1- There then remain only two free parameters which are relevant for 
the dark mater abundance constraint Eq. (1271) : Aj.(= VQ/R') and the bulk mass parameter cn for 
the bulk sterile (right-handed) neutrino. 

For the case of the Planck brane localized Majorana mass, we can, for instance, take (ctv, A^) ~ 
(—0.28, 10 TeV) (corresponding to l/i?' ~ 4 TeV) to satisfy Eq. (I27|l . This in turn gives (m^r. A) = 
(17 keV, 7 X 10^^°) for a warm dark matter candidate, which were obtained by integrating the 5D 
wave functions over the fifth dimension using Eqs. f|T5| [T8|) for the radion decay channel r — > A^A^ 
arising from the term Ar(x)A^(a;)A^(a:;). For the case of the TeV brane localized Majorana mass, 
the choice of (c^r,Ar) ~ (0.63, 1 TeV) (corresponding to 1/i?' ~ 0.4 TeV), for instance, gives us 
{rriN, A) ~ (5 MeV, 4 x 10^^^) so that this case corresponds to a cold dark matter scenario. Note 
we could obtain these small 4D parameters rriisf and A without fine-tuning our free parameters in 
our 5D model thanks to the exponential suppressions coming from the warp factors and small wave 
function overlaps in the extra dimension. 

We also have the constraints from the Lyman-a forest analysis which gives a lower bound for the 
dark matter mass to ensure the enough small scale structures [23j. The Lyman-a forest constrains 
the free-streaming length of the dark matter which in turn gives a constraint on the mass of the dark 
matter particles. The relation between the free-streaming length and the mass however depends on 
the production mechanism affecting the average momentum of the dark matter particles (pn)- For 
a simple estimation of the free-streaming length [1^ |35] 

For the sterile neutrinos produced from the active-sterile neutrino mixing mechanism, (p^) ^ 2.8T 
|46j . while the average momentum for the relativistic fermions in equilibrium is (p^) ~ 3.15T. The 

^For instance, in our simple numerical estimations, we neglected the effects of the the thermal mass corrections 
whose significance would depend on the nature of the electroweak and conformal symmetry breaking phase transitions 



8 



average momentum produced from the thermal radion decay has the average momentum (p^) ^ 
2.45T [19]. Taking account of the additional redshifting of (p^) by the decrease of the effective 
degrees of freedom {g{T ~ TeV)/g{T -C MeV))~^/^, our scenario has {pN)T<^iMeV ~ 0.76T. We 
then find the Lyman-a constraints assuming the the active-sterile neutrino mixing mechanism for 
the sterile neutrino production [23] 

mjv > 10 ~ 28[keV] (29) 

correspond to 0, in our radion decay production scenario, 

"^7v > 3 ~ 8[keV] (30) 

The values realizable in our scenarios, such as those in our concrete examples, hence can be indeed 
consistent with the Lyman-a forest data, and they can also be consistent with another independent 
probe using the QSO gravitational lensing which gives the constraint > 10 keV [47J. 

The other astrophysical constraints, such as the observation of diffuse photon backgrounds 
(211 [25l [271 [291, HH], are related to the neutrino Yukawa coupling and consequently the mixing of 
the active and sterile neutrinos. Our scenarios can satisfy these astrophysical and cosmological 
constraints without affecting the dark matter abundance, simply because the dominant production 
of the sterile neutrinos in our scenarios can come from the radion decay and is unrelated to the 
active-sterile neutrino mixing [211 ISSl [261 1211 HI] (also see Refs. [191 EDj for other production 
mechanisms). In fact, the desirable dark matter abundance is possible even with the vanishing 
neutrino Yukawa coupling for the sterile neutrino in our scenarios . 



6 Discussion and conclusion 

We showed that the sterile neutrino can be an interesting long-lived dark matter candidate in 
a simple warped extra dimension model which in general does not necessarily possess a discrete 
symmetry. 

Although we assumed the (4D) neutrino Yukawa coupling A4/) for the sterile neutrino is negligibly 
small in our discussions, it may be still of possible interest to check how small it has to be. Let us 
for this purpose outline the constraints on the mixing angle sin^ ~ A4/)(if)/mjv where (H) is the 
Higgs vev. 

There arises a constraint from the requirement of the negligible sterile neutrino production by the 
active-sterile neutrino mixings [211 l22l [26] . The fitting function for the sterile neutrino abundance 
produced from the (non-resonant) active-sterile neutrino mixing is0 

""™."'~0-3(^)(l^) (31) 

^°We here, for the number of degrees of freedom, used the Standard Model value {g{T ^TeV)/ g{T ^MeV))^^/'^ ~ 
33-1/3. 

i^The mixing production occurs most efficiently at the temperature around T„iax ~ 133 MeV{mN /IkeVY^"^ 

[211 HI]. 



9 



Another constraint comes from the long enough lifetime [261 ISH EI] by considering the decay into 
an active neutrino and two leptons N u + I + I via Z boson exchange 

which should be at least of order the age of the Universe ~ 4 x 10^^s@ There also exists the 
constraint from the diffuse photon backgrounds [29l [30] due to the radiative decay into an active 
neutrino and a photon N w-f. The emitted monoenergetic photon has a narrow decay line whose 
width is determined by the Doppler broadening for a potential astrophysical probe. For instance, 
for TTiAT < 20 keV, the X-ray observatories such as Chandra and XMM- Newton [23l[25l[27] can give a 
severer constraint than the lifetime constraint. The recent analysis of the XMM- Newton observation 
of the Andromeda galaxy [24j, for example, indicates the mixing angle should be sin^(2^) < 10^^^ 
for one of our examples ~ 17 keV. Even though the more massive cold dark mass cases leading 
to the 7-ray backgrounds still have less precise data for detecting the dark matter decay lines, the 
recent analysis using the high-resolution spectrometer SPI set up on INTEGRAL satellite [58l [59] 
can constrain the dark matter lifetime for 40 keV< rrijsr < 14 MeV and it constrains the mixing 
angle to be less than sin2(26') < lO'^'^ for ~ 5 MeV. ET 

We hence see that our two examples rriN- = 17 keV for a Planck brane localized Majorana mass 
and mjsr = 5 MeV for a TeV brane localized Majorana mass in §3 require the 4D neutrino Yukawa 
couplings A4D which can be obtained after integrating the wave functions over the fifth dimension 
to be respectively A4D < 10~^^ and X^d ^ 10"^'' (we took the Higgs vev (H) = 246 GeV). The free 
parameters in our model to determine X^d are 5D neutrino Yukawa coupling X^d and the bulk mass 
parameters cl for the left-handed neutrino (assuming the other parameter cat relevant for X^d is 
already fixed from the dark matter abundance constraint as in ^5^. It is easy to obtain such a small 
X4D by localizing the left-handed neutrino toward the Planck brane (i.e. by letting cl > 1/2) due 
to the exponential suppressions of the wave function overlaps in the fifth dimension. We however 
found that cl is essentially constrained to be 1/2 < cl < 1 in our simple scenarios to give the 
4D (electron) Yukawa coupling of at least the order O{10^^). This is because the electron and the 
left-handed neutrino are in the same SU(2) doublet and they share the same cl (even though we can 
still freely vary the bulk mass parameter of the right-handed electron). Because of this constraint 
on Cl, to obtain X4^d ^ 10~^^ {X4^d ^ 10"^'') for a warm (cold) dark matter example, we would need 
our free parameter X^d to be X^d ^ 10~^ {X^d ^ 10^^) where the relaxation of tunings compared 
with a 4D model by the order O{10~^^ ~ 10~^) is a result of the exponential suppressions from the 
warp factor and the small wave function overlaps in the extra dimension. 



^^The decay channels of the sterile neutrino via the radion exchange do not affect our discussions because the 
radion interaction vertices are suppressed by S> rriN- 

^•^There are other (less tight) constraints for the heavy sterile neutrinos coming from, for instance, Big Bang 
Nucleosynthesis, supernova and accelerator experiments [SH ISH [551 [Ml and they are satisfied for the small 
mixing angles consistent with the lifetime and photon background constraints. 

^5D} which is an arbitrary free parameter in our model, still however needs to be smaller than the naturally 
expected values of the order unity and its justification goes beyond the scope of the present paper (even though one 
might be tempted to introduce some symmetry for its justification (such as U{1) as discussed for the i^MSM [57]) 
for the vanishing or a small value coming from its breaking). 
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The lightest mass eigenvalue for the light left-handed neutrino in the same generation as that 
of the sterile (dark matter) neutrino is much smaller than O{10^^) eV due to the negligible active- 
sterile neutrino mixing in our model, so that the atmospheric and solar neutrino data [15] should 
be explained by the other two remaining right-handed neutrinos in our simple scenarios. We also 
mention that the desirable baryon asymmetry of the Universe can be induced in our model from 
the decay of the higher KK modes of the right-handed neutrinos [32] . 

While most of the radion phenomenology and in particular the cosmology relevant for the ra- 
dion has been studied for the fields confined on the brane, our study would hopefully open up an 
interesting possibility for the further investigations of both astrophysics and particle physics model 
building in a warped extra dimension with the fields in the bulk. 
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